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The 'Cavity-Mean-Field' approximation developed for the Random Transverse Field Ising Model 
on the Cayley tree [L. loffe and M. Mezard, PRL 105, 037001 (2010)] has been found to reproduce 
the known exact result for the surface magnetization in d = 1 [O. Dimitrova and M. Mezard, J. Stat. 
Mech. (2011) P01020]. In the present paper, we propose to extend these ideas in finite dimensions 
d > 1 via a non-linear transfer approach for the surface magnetization. In the disordered phase, 
the linearization of the transfer equations correspond to the transfer matrix for a Directed Polymer 
, in a random medium of transverse dimension D = d — 1, in agreement with the leading order 

■ perturbative scaling analysis [C. Monthus and T. Garel, arXiv:1110.3145 . We present numerical 

' results of the non-linear transfer approach in dimensions d — 2 and d = 3. In both cases, we find 

^ , that the critical point is governed by Infinite Disorder scaling. In particular exactly at criticality, the 

. one-point surface magnetization scales as Inm^"'-^ ~ —L'^'v, where ujc{d) coincides with the droplet 

' exponent ujdp{D — d — 1) of the corresponding Directed Polymer model, with ajc(d — 2) = 1/3 

and LJc{d = 3) ~ 0.24. The distribution P{v) of the positive random variable v of order 0(1) 
presents a power-law singularity near the origin P{v) oc v" with a{d = 2, 3) > so that all moments 

I of the surface magnetization are governed by the same power-law decay {m^j^^^)^ oc L~^^ with 

5lJ ■ Xs = tJc(l + a) independently of the order k. 
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INTRODUCTION 



The quantum Ising model 



o . . 

^ , where the nearest-neighbor couplings Jij > and the transverse-fields hi > are independent random variables 
drawn with two distributions 'Kcoupiing{J) and T^fieid{h) is the basic model to study quantum phase transitions at 
I zero-temperature in the presence of frozen disorder. In dimension d — 1^ exact results for a large number of observables 
^ , have been obtained by Daniel Fisher [l| via the asymptotically exact Strong Disorder renormalization procedure (for a 
QQ ' review, see Q)- In particular, the transition is governed by an Infinite Disorder fixed point and presents unconventional 
\^ scaling laws with respect to the pure case. In dimension d > 1, the Strong Disorder renormalization procedure can 
still be defined. It cannot be solved analytically, because the topology of the lattice changes upon renormalization, 
\ but it has been studied numerically with the conclusion that the transition is also governed by an Infinite Disorder 
y—{ fixed point in dimensions d = 2, 3, 4 (3l-[l3|. These numerical renormalization results are in agreement with the results 
T-H I of independent quantum Monte-Carlo in d = 2 fl5|. 

^""^ I Even if it is clear that the most natural method to study Infinite Disorder fixed points is the Strong Disorder 
' _ [ ' renormalization approach, it seems useful to determine whether other approaches are able to describe Infinite Disorder 
, scaling. In this paper, we introduce a simple non-linear transfer approximation for the surface magnetization in finite 
dimension d > 1, which is inspired from the 'Cavity- Mean-Field' approximation developed in Refs [T6l - IT8| . and we 
. , study numerically the critical properties of this approximation in dimensions d = 2 and d = 3, 
■ The paper is organized as follows. In Section |lTl we recall briefly the 'Cavity-Mean-Field' approximation developed 
in Refs |16l4l8l | and introduce the non-linear transfer approach for finite dimensions d > 1. Our numerical results in 
dimension d = 2 and d = 3 are presented in sections IIIII and IIVI respectively. Our conclusions are summarized in 
section |Vl 

II. NON-LINEAR TRANSFER APPROACH FOR THE SURFACE MAGNETIZATION 

A. 'Cavity-Mean-Field' approximation on the Cayley tree |l6l 4l8|] 

For the random quantum Ising model model defined on a tree of coordinence {K -\- 1), the following 'Cavity- Mean- 
Field' approximation has been developed [l6l - [l8| : an ancestor i is submitted to the effective single spin Hamiltonian 

Hff = -B^al - h^af (2) 
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where hi is its own random transverse field, and where Bi represents the longitudinal field created by the K children 
j (related to i by the ferromagnetic couplings Jij) within a 'Mean-Field approximation' ( the operator cr| is replaced 
by its expectation value < o-| >) 

K 

The effective Hamiltonian of Eq. [5] is only a two- level system that can be solved immediately : the magnetization of 
the ground state reads 

]3 

=< at >^e//= ' (4) 

V ^, + hi 

Using Eq. |3l one obtains the following non-linear recurrence for the magnetizations rrii (see Eq. 4 of ^Tg*], Eq. 7 of 
[TtI . Eq. 17 of [l^ in the limit of zero temperature /3 = -l-oo) 

. = (5) 



We refer to Refs EMI for more details on this 'Cavity-Mean-Field' approximation and on its properties. As a final 
remark, let us stress that the 'Cavity-Mean-Field' is not exact for the pure model on the Cayley tree (see Fig. 3 of 
Ref [3), but has been argued to become quantitatively correct in the limit of high connectivity if ^ 1 (Tcl - flq . 

In the disordered phase where the magnetizations flows towards zero, the non-linear recurrence of Eq. [5] can be 
linearized to give the following recursion 



1 ^ 

.7 = 1 



J.^jlTlj (6) 



which is equivalent to the problem of a Directed Polymer on the Cayley tree p^| - [l8| . This equivalence can be justified 
directly at the level of lowest-order perturbation theory (i.e. without invoking the 'Cavity- Mean-Field' approximation 
of Eq. [5|), and can be in this way extended to the finite-dimensional case [19|. 



B. 



'Cavity-Mean-Field' approximation in d = 1 [IS 



For K = 1, the Cayley tree of coordinence {K -\- 1) discussed in the previous section becomes a one-dimensional 
chain, and Eq. [5] becomes the one-dimensional non- linear recurrence [l8| 



(Jt^t+im^+i) + hi 



(7) 



Assuming one starts with the boundary condition = 1 at site i = L, one obtains the following explicit expression 
for the surface magnetization nX^^^^ at the site i = jTsI 



surf 



L-1 



hi 



-1/2 



(8) 



As stressed in [18'|, this ex pres sion exactly coincides with the rigorous expression that can be obtained from a free- 
fermion representation (20l. l2l|. and from which many critical exponents can be obtained [2l| - |23| . 

The reason why the 'Cavity-Mean-Field' approximation turns out to become exact for the surface magnetization 
in d = 1 is not clear to us, and seems rather surprising : usually 'mean-field approximation' are exact in sufficiently 
high dimensions or on trees, and are not exact in low dimensions, the 'worst case' being precisely d — I. Here we 
have exactly the opposite conclusion : the 'Cavity-Mean-Field' is not exact for the pure model on the tree (for the 
disordered case, it is not known), but turns out to be exact in d = 1, both for the pure and the disordered case. 
In the absence of any satisfactory explanation for this unusual situation, we tend to think that the exactness of the 
'Cavity-Mean-Field' in c? = 1 is likely to be a 'coincidence' specific to this particular case, from which one cannot draw 
general conclusions for the validity of this approach in higher dimensions or for other quantum disordered models. 
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C. Non-linear transfer approach in finite dimension d > 1 



1. DescnpUon 



periodic b.c. in }^ 




imposed b.c. at x =0 : observables under study: 

m (x -0, y )-l surface magnetizations m(x-L,y) 



FIG. 1: Notations to define tlie non-linear transfer approach in d = 2 : we impose tfie boundary conditions m{x = 0, y) = I on 
tfie left boundary, and we study tlie surface magnetizations m{x — L,y) — 1 on the riglit boundary (see text for more details). 



In finite dimensions d > 1, the authors of Ref [18| have proposed to use the 'Cavity- Mean-Field' approximation 
on a Cayley tree with parameter K = 2d — 1 (to reproduce the connectivity of each spin). In the present paper, we 
propose instead to extend Eqs [5] towards an appropriate non-Unear transfer approach for the surface magnetizations 
of a finite sample of volume L"^. 

For clarity, let us first explain the procedure for the case d = 2. As shown on Fig. [U we consider a lattice containing 
spins : when x is even (x — 0, 2, 4, ..), the coordinate y takes the L integer-values j/ = 1, 2, L ; when x is odd 
{x = 1, 3, ...), the coordinate y takes the L half-integer-values y = 3/2, 5/2, L -f- 1/2. The boundary conditions are 
periodic in y with y + L = y. At x = 0, we impose the boundary condition of unity magnetization 

m{x^O,y)^l (9) 

and we are interested in the L surface magnetizations m{x — L,y) at the opposite boundary x — L. 

For this situation, we propose to use the ideas of Eqs [S] within the following transfer approach. We assume that we 
have already found the surface magnetizations on the column m{x — 1,2/), and we add another column of L sites at 
X. From the Cavity point of view, the new spin at (x, j/) is submitted to its own random transverse field h{x,y) and 
to the longitudinal field (see Eq. |3]) created by its two neighbors on the column x — 1 

B{x,y) = Jf(^.y).(^-i,u+i\xm (x^l,y+l] + Ju,,,-,u.,i\xm [x-l,y-]-] (10) 



'{(x,,),(.-i,,+i)}'" y^^'^y^2)^ ■\{..y).{.-i.y-k)}'" -'^y- 2, 

so that its surface magnetization reads (Eq. 2]) 

/ N B{x,y) 
m(x,y) — — - (11) 

^B^{x,y) + h^ix,y) 

Eqs 1101 and II II define a non-linear transfer procedure that can be iterated from the boundary condition on the column 
a; = of Eq. |9]up to a; = L, where we analyze the statistics of the final surface magnetizations m(x = L,y). It is 
clear that the generalization of this procedure to c? = 3 is straightforward : we add another direction z with periodic 
boundary conditions that plays exactly the same role as y. 
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2. Linearized transfer matrix within the disordered phase 

Within the disordered phase, the surface magnetizations m{x — L,y) are expected to decay typicaUy exponentiahy 
in L, so that one may hnearize the transfer Egs 1101 and 1111 to obtain 



Linearization : mix, y) ~ ——7 r —m [ x — l,y -\ — H r-) r —m { x — l,y (12) 

h{x,y) V ' 27 h{x,y) \ 2J 

This linearized equations can be derived directly within a lowest-order perturbative approach [l9| (i.e. without 
invoking the 'Cavity-Mean- Field' approximation) and corresponds to the transfer matrix satisfied by the partition 
function of a Directed Polymer with D = {d — 1) transverse directions, as discussed in detail in [1^. We refer to 
[TqI for the description of the consequences of this correspondence, and for the analogy with Anderson localization, 
where the droplet exponent of the Directed Polymer also appears in the localized phase [23 - l26| . Here our conclusion 
is that the non-linear transfer approach describes at least correctly the disordered phase, where it coincides with the 
lowest-order perturbative approach jT9| . 



3. Discussion 



Besides its correctness in the disordered phase that we have just discussed, the validity of the non-linear transfer 
exactly at criticality and in the ordered phase has to be studied for the disordered case in d > 1. Since it has been 
found to be exact in d = 1 (see section Hi Bp . one could hope that it is not 'too bad' in d = 2, 3 (even if it is clear 
that this approximation is not valid for the pure model) : we believe that it should capture correctly the nature of 
the transition between 'Infinite-Disorder' or 'Conventional' scaling. In the following, we present our numerical results 
in d = 2 and d = 3 and discuss the scaling properties in the two phases and at criticality. 



III. NUMERICAL RESULTS IN DIMENSION d = 2 



In this section, we present the numerical results obtained with the following sizes L and the corresponding numbers 
rig (L) of disordered samples of volume 

L = lo^ 2.lo^ 3.lo^ 4.lo^ 5.lo^ e.io^, 7.10^, 8.10^ 

n,{L) =2.10^13.10^65.10^37.10^24.10^17.10^13.10^10'' (13) 

For each sample a, we collect the L values of the surface magnetization rrS°'\x = L,i) at the different points 
« = 1, 2, .., L of the surface (see Fig. [ij. Average values and histograms are then based on these L x ns{L) values. 
We have chosen to consider the following log-normal distribution for the random transverse fields hi > 

1 _ (In h-Tnh)^ 

of parameter Inh = and cr = 1, whereas the ferromagnetic couplings Jij are not random but take a single value J 
that will be the control parameter of the quantum transition. 

A. Disordered phase (J < Jc) 

1. Exponential decay of the typical surface magnetization 
In the disordered phase J < Jc, one expects that the typical surface magnetization defined by 



ln(TO^^) = In m^"'^'^ (15) 

decays exponentially with L 

L 



ln(mr)^lnmr^(J<J,)^ --— - (16) 
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FIG. 2: Disordered phase J < Jc in d = 2 (a) Exponential decay of the typical surface magnetization m^"' = e'"'"^ ^ : we 
show the linear decay of In rn^^ as a function of the length L (see Eq. I16p . (b) Log- log plot of the width Ai, of the distribution 
of the logarithm of the surface magnetization as a function of L (here for J = 0.29) : the slope is of order u ~ 0.33 (see Eq. 



where S^typ represents the typical correlation length that diverges at the transition as a power-law 

6,p(J) ^ (J,-J)-'«. (17) 

J->-Jc 

On Fig. [2] (a) we show our numerical results : concerning the exponential decay with L of Eq. [16] for various values 
of J. We find that the corresponding slope l/£,typ{J) vanishes near the critical value Jc ~ 0.335 with the exponent 

Utyp ~ 1 (18) 



2. Growth of the width of the distribution of the logarithm of the surface magnetization 



In the disordered phase J < Jc, one expects that the width of the distribution of the logarithm of the surface 
magnetization defined by 

Ai = ((Inm^'O' - (lnmf^)2) (19) 

grows as a power-law of L 

Ai(J < Jc) ^ (20) 

Our numerical data shown on Fig. [3] (b) correspond to the value 

tj(d = 2)~0.33 (21) 

in agreement with the argument presented in [l^ where uj{d = 2) should coincide with the Directed Polymer droplet 
exponent ujdp{D = d - 1 = 1) = 1/3 



3. Distribution of the logarithm of surface magnetization 

We show on Fig [3] (a) our numerical results concerning histograms of the logarithm of the surface magnetization in 
the disordered phase. Our conclusion is that the surface magnetization follows the scaling 

ln(mf''-^) ~ ln(m^^^') + AiU (22) 

where the behaviors of the typical value ln(m^^'') ~ ~L/£,typ and of the width A^ ^ have been already discussed 
above in Eqs [TOl and respectivelv. On Fig. [3](b), we show that the stable distribution P{u) of the rescaled variable 
u coincides with the GOE Tracy- Widom distribution, as expected from the correspondence with the Directed Polymer 
model in the disordered phase [loj . 
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FIG. 3: Disordered phase in d = 2 (here J = 0.29) : (a) Evolution with L of the probability distribution PL(lnm^'"^''^) of the 
logarithm of surface magnetization : (b) Corresponding fixed distribution of the rescaled variable u = (InmJ"'^-'^ — Inm'j^'') / Al 
in log-scale to show the tails, compared to the exact Tracy- Widom GOE distribution (thick line). 



B. Ordered phase 

1. Behavior of the typical surface magnetization in the ordered phase J > Jc 




0.33 0.34 0.35 0.36 0.37 0.38 



FIG. 4: Ordered Phase J > Jc m d — 2 : (a) Behavior of the asymptotic typical surface magnetization m}^' as a function of 
the ferromagnetic coupling J : our numerical data are compatible with an essential singularity (Eq. I24p of exponent k ~ 0.5 
(b) the probability distribution Pt(lnm^'"^''^) (here for J — 0.38) of the logarithm of surface magnetization remains fixed and 
attached at the origin (as a consequence of the bound m^"'^^ < 1). 



In the ordered phase, the typical surface magnetization remains finite in the limit where the number of generations 
L diverges 



lnm^^P(J > Jc) = lnm''^''^{J > Jc) ~ lnTOoo(J > Jc) > -oo (23) 

L— >oo 

and one expects an essential singularity behavior 

\nm'y^P{J > Jc) oc -(J-Jc)-^ (24) 

Our data shown on Fig. S] (a) can be fitted with the value 

K(d = 2)~0.5 (25) 
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that can be related to other exponents via finite-size scahng (see below around Eq. l34l) 



2. Distribution of the logarithm of the surface magnetization 

In the ordered phase, the probability distribution Pi(lnm^"'^''^) of the logarithm of surface magnetization remains 
fixed as L varies, and terminates discontinuously at the origin, as a consequence of the bound mf^j^^^ < 1 corresponding 
to Inmf''^ < (see Fig. H (b)) 

C. Critical point 

1. Behavior of the typical surface magnetization rn^^ and of the width Ai 
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FIG. 5: Critical point (here Jc = 0.335) : (a) Log-log plot of the logarithm of the typical surface magnetization m^"' and of the 
width A_L : both slopes are of order Wc — 0.33 (see Egs 1261 and I27|) (b) Finite-size scaling of the typical surface magnetization 
according to Eq. [29]with Jc = 0.335, ujc = 0.33 and Uav ~ 1-5. 

Exactly at criticality, one expects that the typical surface magnetization follows an activated behavior of exponent 
cjc < 1 (compare with Eq. \W\ in the disordered phase) 



ln(m^^''(J = Jc))=lnmf''-^(J = Jc) c 

and that the width defined in Eq. [T5] is also governed by the same exponent 

Ai(J = J,) ~ L^^ 

L-^oo 



(26) 



(27) 



Our numerical data at Jc ~ 0.335 shown on Fig. [5] (a) are compatible with these behaviors with the value 

ujc ~ 0.33 (28) 

i.e. ujc coincides with the fluctuation exponent u measured in the disordered phase (see Eq. I2ip 

This last property implies that the finite-size scaling for the typical surface magnetization to^*^ involves some 
correlation length exponent Vav different from vtyp 

Inm^^P(J) = In mf^(J)~-L"=G = ( J - J^L^/''"') (29) 

The matching with the behavior of Eq. [55] in the disordered phase implies that 

G{x) cx i-xY""' (30) 
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and that Vav reads 

^a. = ^ (31) 

1 — w 

This relation can be understood within a rare events analysis for the averaged correlation in the disordered phase 
The values ^typ — 1 and w — 1/3 yield 

VaAd^2)c^^ (32) 

The matching of the finite-size scaling form of Eq. with the essential singularity of Eq. [M] in the ordered phase 
implies that 

G{x) cx \ (33) 

with 

K = UJc^av = l^typ-r—— (34) 

<d = 2) ~ i (35) 

in agreement with the estimate of Eq. [511 

As shown on Fig[5] (b), our numerical data collapse well with the finite-size scaling form of Eq. [29] with Vav — 1-5. 



The values vtyp = 1 and w = 1/3 yield 



2. Distribution of the logarithm of surface magnetization 




FIG. 6: Critical point in d = 2 (here Jc = 0.335) 
logarithm of surface magnetization. Inset : Corresponding fixed distribution of the rescaled variable v — {Iw: 



Evolution with L of the probability distribution PL(him^j^^^) of the 



surf 
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At criticality, the rescaled variable 



1 

in m 



surf 



In 771 



surf 



In m^^^ 



(36) 



remains a positive random variable of order 0(1) as L — > +oo. Our numerical measure of its probability distribution 
P{v) shown on Fig. |6]is compatible with a power-law singularity near the origin 



P{v) 



(37) 



with an exponent of order a > 2 that we do not measure precisely. Note that this is different from the case d = 1 
where P{v = 0) is finite (a = 0). We have not been able to find a physical argument to predict the value of a in d = 2. 
This exponent a will directly influence the scaling of the moments of the surface magnetization, as we now discuss. 



9 



3. Moments of the surface magnetization 

In contrast to the activated behavior of the typical surface magnetization to^^ of Eq. I26[ the moments of the 
surface magnetization are expected to fohow a power-law, as a consequence of the following rare events analysis : the 
surface magnetization of Eq. |3S]will be of order 0(1) if the random variable v happens to be smaller than l/L'^". 
Taking into account the behavior of Eq. [37l this will happen with probability 

Prob{m''^''f dvP{v) - / dvv" cx (38) 







and all moments will be governed by this power-law 



(mf ~ Pro6(mf "-^ = 1) oc with = + a) (39) 

independently of the order k. 

Our numerical data for the three first moments k — 1,2,3 and various sizes are compatible with Eq. 1391 with an 
exponent 

a;,(d = 2)~1.2 (40) 

The relation of Eq. 131] then corresponds to 

a(d=2)~2.6 (41) 
In the ordered phase, our numerical data are compatible with the power-law 



mf^ cx ( J - Jcf' (42) 

with 

j3s{d^2)^XsVav-^l-^ (43) 



IV. NUMERICAL RESULTS IN DIMENSION d = 3 



In this section, we present the numerical results obtained with the following sizes L and the corresponding numbers 
ns{L) of disordered samples of volume 

L = 10^ 2.10^ 3.10^ 4.10^ 5.10^ e.lO^ 7.10^ 8.10^ 
UsiL) = 27.10^7.10^ 3. 10^ 16. 10^10^ 7.10^,5.102,4.102 (44) 

For each sample a, we collect the values of the surface magnetization at the different points of the surface. Average 
values and histograms are then based on these x ns{L) values. We consider again the disorder distribution of Eq. 
[13] and take J as the control parameter of the transition. 



A. Disordered phase 

Our data follow the scaling of Eq. |22l with the following properties : 

(i) the scaling of the typical surface magnetization is given by Eq 1161 ^-nd the typical correlation length exponent 
of Eq. [TT] seems again very close to unity 

Vtyp ~ 1 (45) 

(ii) the width of Eq. [12] grows as the power-law of Eq. [^D] with the exponent 

uj{d = 3)~ 0.24 (46) 

that coincides with the numerical values of the droplet exponent of the Directed Polymer model with D = d — 1 = 2 
transverse dimensions (3ll - [37j . in agreement with the argument presented in [l9j . 

(iii) As L grows, the evolution of the probability distribution PL(lnTO^"'^'^) is shown on Fig. [7] (a). The corresponding 
fixed distribution of the rescaled random variable u = (Inm^"'^^ — lnm^'')/Ai is shown on Fig [7] (b). 
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FIG. 7: Disordered phase in d = 3 (here J = 0.11) : (a) Evolution with L of the probability distribution PL(lnm^'"^''^) of the 
logarithm of surface magnetization : (b) Corresponding fixed distribution of the rescaled variable u = (InmJ"'^-'^ — Inrn^^) / 



in log-scale to show the tails. 



B. Critical point 




FIG. 8: Typical surface magnetization in the critical region in d = 3 (a) Log-log plot of the logarithm of the typical surface 
magnetization rn^^ and of the width A_t at criticality Jc — 0.1528 : both slopes are of order cuc — 0.24 (see Eqs 1261 and 1271 ) 
(b) Finite-size scaling of the typical surface magnetization in d = 3 according to Eq. [29]with i^av = 1.32. 



At criticality Jc — 0.1528, we find that the exponent Wc of Eqs [55] and [57] coincides with the value of ui of Eq. [JS] 
concerning the disordered phase (see Fig. [5] (a)) 

uJc{d = 3) ~ 0.24 (47) 

We show on Fig. [H](b) the finite-size scaling analysis of Eq. [2S]for the logarithm of the typical surface magnetization 
with the averaged correlation length exponent Vav = 1/(1 — w) ~ 1.32. 

We show on Fig. [5] our numerical data for the probability distribution of the surface magnetization at criticality : 
the fixed point distribution P{v) of the rescaled variable v of Eq. [35] displays a power-law singularity near the origin 
(Eq 1571) , that will determine the scaling of all moments of the surface magnetization according to Eq. [3S] Our 
numerical data for the moments are compatible with Eq. [39] with an exponent 

a;s(d = 3) ~ 1.34 (48) 

that would correspond to 



a((i = 3) = — - 1 ~ 4.5 (49) 

OJc 
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FIG. 9: Critical point in d = 3 (here Jc — 0.1528) : Evolution with L of the probability distribution PL(lnm^"^ ) of the 
logarithm of surface magnetization. Inset : Corresponding fixed distribution of the rescaled variable v — (In mj^""^^)/ Inm^^''. 

and to the exponent (Eq. 

I3s{d = 3) = xsiyav^ 1-76 (50) 
V. CONCLUSION 

Since the 'Cavity-Mean-Field' approximation developed for the Random Transverse Field Ising Model on the Cayley 
tree has been found to reproduce the known exact result for the surface magnetization in d = 1 [isj . we 

have proposed to extend these ideas in finite dimensions d > 1 via a non-linear transfer approach for the surface 
magnetization. In the disordered phase, the linearization (Eg I12p of the transfer equations correspond to the transfer 
matrix for a Directed Polymer in a random medium of transverse dimension D = d — I, in agreement with the leading 
order perturbative scaling analysis [l9j . 

We have presented numerical results of this non-linear transfer approach in dimensions d = 2 and d = 3, where 
large system sizes can be easily studied. In both cases, we have found that the critical point is governed by Infinite 
Disorder scaling. In particular exactly at criticality, the one-point surface magnetization scales as Inm^"'^-'^ ~ —L'^''v, 
where Ludd) coincides with the droplet exponent ujdp{D = d — 1) of the corresponding Directed Polymer model, 
with ujc{d = 2) = 1/3 and uJc{d = 3) ~ 0.24. The distribution P{v) of the positive random variable v of order 0(1) 
presents a power-law singularity near the origin P{v) oc w° with a{d = 2,3) > so that all moments of the surface 

magnetization are governed by the same power-law decay {rn^j^^^)^ cx L^^= with Xg = Wc(l + a) independently of the 
order k. Our conclusion is thus that this non-linear transfer approach is able to lead to Infinite Disorder scaling, that 
had been found previously via Monte-Carlo in d = 2 [l3,[l^ and via Strong Disorder RG in d = 2, 3, 4 Exactly 
at criticality, the presence of activated scaling Inm^"'''^ ~ —L^'v means that the linearization of Eq. [T2]is typically 
still valid also at criticality (and not only in the disordered phase), so that the identity ujc = ujdp{D = d — 1) can be 
understood. The rare cases where this linearization is not valid at criticality is when the positive random variable v 
happens to be smaller than 1 /L'^" . Our conclusion is thus the following : 

(i) in the disordered phase and for 'typical' situations exactly at criticality, the linearization of Eq. [T2]is valid and 
coincides with the leading order perturbative scaling analysis [l^ : it is thus likely to give exact values for critical 
exponents, in particular for the exponent uJc of activated scaling. 

(ii) in the ordered phase and for 'rare' situations at criticality, the non-linear terms of the transfer approach plays 
an important role. Since they come from an uncontrolled approximation, the critical exponents like Ps and Xs that 
are determined by these non-linear contributions could be different from the exact ones. To judge the accuracy of this 
approximation, it would be very helpful to compare with other approaches like Quantum Monte-Carlo and Strong 
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Disorder RG (but up to now, these other approaches have not studied surface properties). 



[1] D. S. Fisher, Phys. Rev. Lett. 69, 534 (1992); Phys. Rev. B 51, 6411 (1995). 
[2] F. Igloi and C. Monthus, Phys. Rep. 412, 277 (2005). 
[3] D. S. Fisher, Physica A 263, 222 (1999). 

[4] O. Motrunich, S.-C. Mau, D. A. Huse, and D. S. Fisher, Phys. Rev. B 61, 1160 (2000). 

[5] Y.-C. Lin, N. Kawashima, F. Igloi, and H. Rieger, Prog. Theor. Phys. 138, 479 (2000). 

[6] D. Karevski, YC Lin, H. Rieger, N. Kawashima and F. Igloi, Eur. Phys. J. B 20, 267 (2001). 

[7] Y.-C. Lin, F. Igloi, and H. Rieger, Phys. Rev. Lett. 99, 147202 (2007). 

[8] R. Yu, H. Saleur, and S. Haas, Phys. Rev. B 77, 140402 (2008). 

[9] I. A. Kovacs and F. Igloi, Phys. Rev. B 80, 214416 (2009). 
[10] I. A. Kovacs and F. Igloi, Phys. Rev. B 82, 054437 (2010). 
[11] I. A. Kovacs and F. Igloi, Phys. Rev. B 83, 174207 (2011). 
[12] I. A. Kovacs and F. Igloi, arXiv:1108.3942 
[13] I. A. Kovacs and F. Igloi, arXiv: 1109.4267 

[14] C. Rich, A. P. Young, H. Rieger, and N. Kawashima, Phys. Rev. Lett. 81, 5916 (1998). 

[15] H. Rieger and N. Kawashima, Eur. Phys. J B9, 233 (1999). 

[16] L. B. loffe and M. Mezard, Phys. Rev. Lett. 105, 037001 (2010) 

[17] M. V. Feigelman, L. B. loffe, and M. Mezard, Phys. Rev. B 82, 184534 (2010). 

[18] O. Dimitrova and M. Mezard, J. Stat. Mech. (2011) P01020. 

[19] C. Monthus and T. Garel, arXiv:1110.3145 

[20] I. Peschel, Phys. Rev. B 30, 6783 (1984). 

[21] F. Igloi and H. Rieger, Phys. Rev. B 57, 11404 (1998). 

[22] A. Dhar and A.P. Young, Phys. Rev. B 68, 134441 (2003). 

[23] C. Monthus, Phys. Rev. 69, 054431 (2004). 

[24] V.L. Nguyen, B.Z. Spivak and B.I. Shklovskii, JETP Lett. 41, 43 (1985); V.L. Nguyen, B.Z. Spivak and B.I. Shklovsku, 

Sov. Phys. JETP 62„ 1021 (1985). 
[25] E. Medina, M. Kardar, Y. Shaph and X.R. Wang, Phys. Rev. Lett. 62, 941 (1989); E. Medina and M. Kardar, Phys. Rev. 

B 46, 9984 (1992). 

[26] J. Prior, A.M. Somoza and M. Ortuno, Phys. Rev. B 72, 024206 (2005); A.M. Somoza, J. Prior and M. Ortuno, Phys. 

Rev. B 73, 184201 (2006); A.M. Somoza, M. Ortuno and J. Prior, Phys. Rev. Lett. 99, 116602 (2007). 
[27] D. A. Huse, C. L. Henley, and D. S. Fisher, Phys. Rev. Lett. 55, 2924 (1985). 
[28] M. Kardar, Nucl. Phys. B290, 582 (1987). 
[29] K. Johansson, Comm. Math. Phys. 209, 437 (2000). 

[30] M. Prahofer and H. Spohn, Physica A279, 342 (2000); M. Prahofer and H. Spohn, Phys. Rev. Lett. 84, 4882 (2000); M. 

Prahofer and H. Spohn, J. Stat. Phys. 108, 1071 (2002); M. Prahofer and H. Spohn, J. Stat. Phys. 115, 255 (2002). 
[31] L.H. Tang, B.M. Forrest and D.E. Wolf, Phys. Rev. A 45 (1992) 7162. 
[32] T. Ala-Nissila, T. Hjeff, J.M. Kosterlitz and V. Venalainen, J. Stat. Phys. 72 (1993) 207. 
[33] E. Perlsman and M. Schwartz, Physica A 234, 523 (1996). 

[34] T. Ala-Nissila, Phys. Rev. Lett. 80 (1998) 887 ; J.M. Kim, Phys. Rev. Lett. 80 (1998) 888. 

[35] E. Marinari, A. Pagnani and G. Parisi, J Phys. A33, 8181 (2000); E. Marinari, A. Pagnani and G. Parisi and Z. Racz, 

Phys. Rev. E65, 026136 (2002). 
[36] C. Monthus and T. Garel, Phy s. Rev. E 73 , 056 106 (2006); C. Monthus and T. Garel, Phys. Rev. E 74, 051109 (2006). 
[37] M. Schwartz and E. Perlsman, larXiv: 1108.4604] 



